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In the determination of ttj from tau decays several different moments of the hadronic spectral 
I""! functions have been used. In a recent work, we performed an analysis of their perturbative be- 

^ We showed that the various moments can be divided in a small number of classes. We con- 

D eluded that some of the moments commonly employed in a, extractions should be avoided due to 

their bad perturbative behaviour. Furthermore, for the moments that have a good perturbative be- 
haviour, and under reasonable assumptions for the higher-order behaviour of the Adler function, 
J> fixed-order perturbation theory (FOPT) provides the superior framework for the renormalization 

group improvement. Here we discuss an extension of this analysis where we consider the pertur- 
<^ bative series for values of the hadronic invariant mass squared SQ<rn^. Our conclusions are not 

altered within a reasonable window. 
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1. Introduction 

In the last 20 years, hadronic tau decays have been an important source of empirical informa- 
tion on fundamental parameters of QCD. Notably, the strong coupling, a^, can be extracted with 
a good precision at relatively low energies, close to the edge of the validity of perturbative QCD. 
After the measurement of the spectral functions at LEP, other parameters such as the strange quark 
mass, the CKM matrix element Vi„, as well as non-perturbative condensates could be extracted 
(see e.g. [1]). The extraction of these parameters relies on sum rules. Quark-hadron duality and the 
optical theorem allow to express the decay rate as weighted integrals of the vector and axial-vector 
spectral functions running over the hadronic invariant mass squared s from threshold up to m|. 

Using the analytic properties of the quark-antiquark correlators, the theoretical counter-part 
of the experimental quantities are expressed as contour-integrals in the complex energy plane with 
fixed l^l = m\. However, in the theoretical description of T decays two main obstacles remain. 
The first is referred to as duality violations (DVs). They are related to the break-down of local 
quark-hadron duality in the vicinity of the Minkowski axis (where resonance effects become im- 
portant). In the past, they have been neglected due to a fortuitous kinematical suppression of the 
problematic region in the contour integration. Recently, thanks to the progress in modelling DVs 
realistically [3, 4], they have been included self-consistently in a full-fledged analysis [5, 6]. The 
second important obstacle is the prescription for the renormalization group (RG) improvement of 
the perturbative series. The most widely employed prescriptions are fixed-order perturbation theory 
(FOPT) [7, 8] and contour improved perturbation theory (CIPT) [9, 10]. When used in practice, 
they lead to different a, results. With the recently computed correction [1 1], the difference be- 
came even more pronounced. Several works have dealt with this discrepancy [8, 12, 13, 14, 15, 16] 
in the light of the term. The conclusions in favour of FOPT or CIPT (or a third prescription) 
are based on (implicit or explicit) assumptions on the yet unknown higher order q corrections. In 
this context, the goal of Ref . [8] was to construct a plausible model for the higher-order connections 
of the Adler function from the leading renormalon singularities of its Borel transform, using only 
general RG arguments to describe the structure of the singularities in the Borel plane. After match- 
ing the model to the known coefficients in QCD, the main conclusion of Ref. [8] was that FOPT is 
to be preferred over CIPT, since FOPT provides a closer approach to the Borel resummed results 
— in the spirit of an asymptotic series. 

This conclusion was based solely on the analysis of the weight Wx{x), obtained from the kine- 
matics of the decay. This is not entirely satisfactory since realistic determinations of q employ 
(and often require) several different weight functions Wi{x). In fact, any analytical Wi{x) gives 
rise to a valid sum-rule that emphasises a given part of the spectral functions, as well as different 
contributions in the theoretical description. In the literature, several weight-functions have been 
employed and yet little attention has been paid to the moment dependence of the convergence 
properties of the perturbative series. We have addressed this question in Ref. [17] and pursued 
the FOPT/CIPT comparison for several weight functions. We showed that the different moments 
employed in the literature can be divided in a small number of categories. The characteristics of 
their perturbative series could be linked to generic features of the moment weight function and the 
dominant renormalon singularities of the Adler function. We concluded that some of the moments 
currently employed in some q extractions should be avoided due to the poor convergence of their 
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perturbative expansions. Additionally, for all moments that display good perturbative behaviour 
— and under reasonable assumptions for the higher-order corrections — FOPT provided the best 
framework to the RG improvement. 

A point that was not discussed in [17] is the stability of the conclusions with respect to 
variations (^o being the upper limit of integration in the sum-rule, see Eq. (2.1)). The relevance of 
this issue lies in the fact that several Us analyses use sum-rules where the data are integrated up to 
sa < m^. Here we show that the conclusions of [17] remain valid when is varied away from 



2. Theoretical framework, model, and results 

We work with generalized sum-rules, where the weight function in the integrals can be any 
analytical function Wi{x) and the upper limit of integration is taken to be any point sq < m^. The 
experimental side of the sum-rules are then written as integrals over the spectral functions as 



/ \ 2 / 

ds I s \ I ^ s 



ImUilr'Hs) ?^Imn!l(.) 



(2.1) 

The two point functions are defined as n(^/^(p) = iJdxe'P'{a\ T{J^^^{x)j;^^^{Oy}\Q.) and they 
assume the usual decomposition into longitudinal and transversal components. The V and A cur- 
rents are given by Jy^^{x) = {Hy^ {Y5)d){x). 

The theoretical counter-part of Eq. (2.1) is obtained from the counter-clock wise contour inte- 
gration of the correlators with \s\ = sq. The contributions on the theory side can be organized as 



D>2 



(2.2) 



The perturbative contribution is contained in 5^*^^ and 5,,- , of which , contains the loop correc- 
tions. In the chiral limit they are the same for V and A correlators, and correspond to the perturba- 
tive series of 11^^^°^ {s) . The quark-mass corrections, as well as contributions from operators with 

D > 2 in the OPE, are encoded in the terms S^^y/j^', DV contributions are represented by 5^^/^- 

Our focus is on the behaviour of the perturbative correction and it is convenient to write it in 
terms of the RG invariant Adler function, whose expansion in can be written as 

z><.«.M . -.^n"«.(,) = ^.i'.ip^.. (.og^)" , (2.3) 

where = a{pL)/n. RG invariance implies that only the coefficients c„j are independent. The 
other Cn^k can be expressed in terms of c„.i and jS -function coefficients. The perturbative contribu- 
tion to the theory side of the sum-rules is then 

= £ £fcc,a ^. f -W,-(x)log^-if ^"ja;^, (2.4) 

1^1=1 

with X = s/sq and Wi{x) = 2 dzwi{z)- Due to the RG invariance of d('+°) (s) one has the freedom 
of setting the scale jj.. The FOPT prescription corresponds to the choice /x^ = sq. In this case. 
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the coupling a{sQ) is taken out-side the integrals and one is left with the integration of powers of 
log(— x). The CIPT choice correspond to /i^ = —sqx, which resums the logarithms but, in turn, the 
integrals are done (numerically) over the running coupling a{—sax). 

In order to compare FOPT and CIPT as well as understand the perturbative behaviour of 
spectral function moments, one must have an ansatz for the unknown higher-order Adler function 
coefficients c„^i . Here we follow the method introduced in Ref. [8] which makes use of the available 
knowledge of the renormalon structure of the Borel transformed Adler function. The idea is to 
construct a realistic model for the Borel transform using the leading singularities. We work with 
the function D{s) and its Borel transform, B[D\{t), defined as 

-^d['+'\s) ^ l+D{s) ^ l + £r„a,(V^r+^ B[D]{t) ^ £r„-. (2.5) 
The original series can be understood as an asymptotic expansion of the inverse of B[D] (t), 

noo 

D{a)= dte-'^^BlD^t), (2.6) 
Jo 

when the integral exists. Singularities of B[D]{t) on the positive real axis (infra red (IR) renor- 
malons) give rise to fixed-sign asymptotic series and obstruct the Borel summation, Eq. (2.6). This 
introduces an ambiguity in the integral that is expected to be cancelled against exponentially small 
terms in a,, or power corrections(due to the logarithmic running of the coupling). Singularities on 
the negative real axis (ultra violet renormalons (UV)) give rise to sign- alternating series. 

General RG arguments and the structure of the OPE allow one to determine the position and 
strength of the renormalon singularities in the t plane, though not their residues [18]. The fixed-sign 
nature of the exactly known coefficients of the Adler function suggest that at low and intermediate 
orders the series is dominated by IR singularities. The reference model (RM) of [8] contains the 
first two IR and the leading UV singularities. The Borel transform is given by 

B[D]iu) = B[Dn(") +B[Df]{u) +B[Df]{u) + (2-V) 

The structure of the branch-cut singularities can be found in [8] . The residues and the coefficients 
d^'^ are fixed by matching to the exactly known ci j to C4 i (augmented by an estimate for cj 1). 

Within this model, the conclusion of Ref. [8] in favour of FOPT has been corroborated and 
extended in our recent work [17]. All moments that display a good perturbative behaviour favour 
the FOPT prescription within the RM. This conclusion can be traced back to the contribution of the 
leading IR singularity, related to the D = 4 corrections in the OPE. If this singularity is arbitrarily 
suppressed, one generates a model — less realistic, in our opinion — in which CIPT is the preferred 
prescription. To realize this scenario in practice, and assess possible model dependencies in our 
conclusions, we introduced the following alternative model (AM) where the leading singularity is 
absent whereas the sub-leading one at m = 4 is explicitly taken into account: 

B[D]{u) = B[D^^]{u)+B[Df]{u)+B[Df]{u) + d^° + d\^u. (2.8) 

Within the AM, moments with good perturbative behaviour favour CIPT. 
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Figure 1: Reference model. (.so) order by order in normalised to the Borel sum for FOPT (left) and 
CIPT (right) with three values of 1-5 GeV^, 2.5 GeV^, and m^. Bands give the Borel ambiguities. 

The models represent two quite different situations regarding the interplay of the Adler func- 
tion coefficients and the running coupling effects. In the RM, there are cancellations between the 
contribution from the high-order coefficients c„ i and the running coupling effects, at a given order 
in ttv In this case FOPT is superior since it treats these contributions on an equal footing, while 
CIPT misses the cancellations due to the resummation of the running effects to all orders. On the 
other hand, the AM represents a situation where the running effects are dominant and should be 
resummed. In this case, the high-order coefficients can be neglected and CIPT is a better prescrip- 
tion. Since there is no known mechanism that would naturally suppress the leading IR singularity 
in QCD, we believe the scenario of Eq, (2.7) to be more realistic. 

Using these two models, we compared in Ref. [17] the perturbative series in FOPT and CIPT 
generated from 17 polynomial weight-functions Wi{x). We showed that they can be divided into a 
small number of categories regarding the behaviour of their perturbative series. Generic features 
of the functions Wi{x) (such as starting or not with the unity), together with the assumptions upon 
the Adler function, suffice to determine whether they are suitable for as extractions and whether 
FOPT or CIPT is more suitable for the RG improvement. We showed that some of the weight 
functions used in the literature, e.g. polynomials containing solely powers of x' with / > 2, should 
be avoided due to their bad perturbative behaviour. We also provided further arguments that support 
the plausibility of the RM of [8] and concluded that for well-behaved moments FOPT is preferred. 

An aspect that was not considered in [17] was the sq dependence of these conclusions. This is 
important because sum-rules with different values of < 'w^ are used in extractions of [5, 6, 19]. 
Here we show explicit results for the FOPT/CIPT comparison for two moments within the two 
models given in Eqs. (2.7) and (2.8) and considering three values of 1-5 GeV^, 2.5 GeV^, and 
ni^. The interval [1.5 GeV^: m^] spans the values used in the fits of [5, 6]. Since we intent to 
compare the perturbative series at different values of a normalisation procedure is in order. For 
better comparison, we normalise the series generated for each value of by its corresponding 
Borel sum, Eq. (2.6). Hence, in the plots, meaningful series should be asymptotic to the unity. 

We start by considering the case of moments that have good perturbative behaviour for 5^0 = ni^- 
As a representative we choose to use the kinematic moment w^- In Fig. 1, we consider the FOPT 
and CIPT series within the RM. On the left-hand side. Fig. 1(a), one observes that the normalised 
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Figure 2: Reference model. Svv,, (sq) order by order in a, normalised to the Borel sum for FOPT (left) and 
CIPT (right) with three values of sq^ 1-5 GeV^, 2.5 GeV^, and m^. Bands give the Borel ambiguities. 

FOPT series still behaves as a good asymptotic series even for significantly smaller than m^. 
As expected, for lower ^o, the larger values of a, amplify the divergent behaviour above the 8th 
order. Nevertheless, the first few terms of the series approach the Borel resummed value. Note 
also that the Borel sum has a larger ambiguity for smaller due to larger as- On the right-hand 
side, in Fig. 1(b), one sees that the poor performance of CIPT is amplified by the larger values of 
the coupling at lower sq. That is, within the RM, CIPT is not a good approximation to the Borel 
resummed values, and even less so for smaller sq. 

We now turn to a moment with bad perturbative behaviour: w\i{x) = (1 — ;c)^.x^(l +2.x) (to 
employ the notation of [17]). In Ref. [17] we showed that moments starting with powers of x 
(that do not contain the unity) tend to have bad perturbative behaviour and are largely dominated 
by power corrections. In Fig. 2 we address the s^^ dependence of this conclusion. On the left. 
Fig. 2(a) shows that for higher values of sq FOPT can approach the Borel result only at high orders 
(not available exactly). At low sq the series displays a wild behaviour and cannot be considered a 
good approximation to the Borel sum. In CIPT, Fig. 2(b), the bad behaviour already observed for 
sq = m\ is amplified at lower ^o- The series are erratic and cannot be consider suitable asymptotic 
approximations to the Borel sum. Note that this moment, despite of its bad perturbative behaviour, 
enters several determinations of from T decays (e.g. Refs. [20, 21]). 

We can perform the same analysis in the alternative model, Eq. (2.8), which receives no contri- 
bution from the leading IR singularity. In Fig. 3, we show the series normalised to their respective 
Borel resummed values within the AM for Wxix). In this model, CIPT provides the better frame- 
work also for lower values of sq. as shown in Fig. 3(b). The series remains very stable for = 2.5 
GeV^ and still approaches the Borel result well. For = 1-5 GeV^ CIPT can still be considered a 
good approximation taking into account the amplified Borel ambiguity. The oscillations of FOPT, 
already present at sq = m\, are much amplified for lower sq (see Fig. 3(a)). Within the AM, the 
FOPT series are not a good approximation to the Borel resummed values. 

To conclude we examine the case of wn in the context of the AM. The results ai^e shown in 
Fig. 4. The bad perturbative behaviour of FOPT and CIPT remains for all values of sq. This is an 
indication of the model independence of the conclusion that wn (and a number of other moments 
also discussed in Ref. [17]) should be avoided in determinations of a^. 
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Figure 3: Alternative model. 5wJ (so) order by order in a, normalised to the Borel sum for FOPT (left) and 
CIPT (right) with three values of sq. 1.5 GeV^, 2.5 GeV^, and m^. Bands give the Borel ambiguities. 
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(a) wn(x) = (l-x)V(l+2x), FOPT. (b) wi7(j) = (1 -x) V(l +2ji;), CIPT. 

Figure 4: Alternative model. S,!*)' {sq) order by order in normalised to the Borel sum for FOPT (left) and 
CIPT (right) with three values of sq: 1.5 GeV^, 2.5 GeV^, and m^. Bands give the Borel ambiguities. 



3. Conclusions 

Recently, we have analysed the perturbative behaviour of several moments often employed 
in analyses of a, from T decays under different assumptions for the large-order behaviour of the 
Adler function [17]. We have shown that some of these moments should be avoided due to their 
bad perturbative behaviour. Furthermore, under reasonable assumptions for the Borel transformed 
Adler function, we showed that FOPT provides the preferred framework for the RG improvement 
of moments that display good perturbative behaviour. 

Here we showed, for the first time, that these conclusions are still valid if one considers the 
perturbative series generated by FOPT and CIPT for sq < m^. This is a relevant question, since in 
as extractions one often considers sum-rules with sq < m^. We have shown explicitly the results for 
two representative moments previously investigated in Ref. [17] for sq = m^. The dependence 
analysis was also earned out for the remaining moments studied in [17] with similar conclusions; 
they are not shown here for the sake of brevity. 
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